Selected Nonorthogonal Configuration Interaction with Compressed Single and Double Excitations.

Addressing both dynamic and static correlations accurately is a primary goal in electronic structure theory. Nonorthogonal configuration interaction (NOCI) is a versatile tool for treating static correlation, offering chemical insights by combining diverse reference states. Nevertheless, achieving quantitative accuracy requires the inclusion of the missing dynamic correlation. This work introduces a framework for compressing orthogonal single and double excitations into a NOCI of a much smaller dimension. This compression is repeated with each Slater determinant in a reference NOCI, resulting in another NOCI that includes all of its single and double excitations (NOCISD), effectively recovering the missing dynamic correlations from the reference. This compressed NOCISD is further refined through a selection process using metric and energy tests (SNOCISD). We validate the effectiveness of SNOCISD through its application to the dissociation of the nitrogen molecule and the hole-doped two-dimensional Hubbard model at various interaction strengths.

Hartree-Fock-Bogoliubov theory for number-parity-violating fermionic Hamiltonians.

It is usually asserted that physical Hamiltonians for fermions must contain an even number of fermion operators. This is indeed true in electronic structure theory. However, when the Jordan-Wigner (JW) transformation is used to map physical spin Hamiltonians to Hamiltonians of spinless fermions, terms that contain an odd number of fermion operators may appear. The resulting fermionic Hamiltonian thus does not have number parity symmetry and requires wave functions that do not have this symmetry either. In this work, we discuss the extension of standard Hartree-Fock-Bogoliubov (HFB) theory to the number-parity-nonconserving case. These ideas had appeared in the literature before but, perhaps for lack of practical applications, had, to the best of our knowledge, never been employed. We here present a useful application for this more general HFB theory based on coherent states of the SO(2M + 1) Lie group, where M is the number of orbitals. We also show how using these unusual mean-field states can provide significant improvements when studying the JW transformation of chemically relevant spin Hamiltonians.

Correlated pair ansatz with a binary tree structure.

We develop an efficient algorithm to implement the recently introduced binary tree state (BTS) ansatz on a classical computer. BTS allows a simple approximation to permanents arising from the computationally intractable antisymmetric product of interacting geminals and respects size-consistency. We show how to compute BTS overlap and reduced density matrices efficiently. We also explore two routes for developing correlated BTS approaches: Jastrow coupled cluster on BTS and linear combinations of BT states. The resulting methods show great promise in benchmark applications to the reduced Bardeen-Cooper-Schrieffer Hamiltonian and the one-dimensional XXZ Heisenberg Hamiltonian.

Symmetry-projected cluster mean-field theory applied to spin systems.

We introduce Sz spin-projection based on cluster mean-field theory and apply it to the ground state of strongly correlated spin systems. In cluster mean-fields, the ground state wavefunction is written as a factorized tensor product of optimized cluster states. In previous work, we have focused on unrestricted cluster mean-field, where each cluster is Sz symmetry adapted. We here remove this restriction by introducing a generalized cluster mean-field (GcMF) theory, where each cluster is allowed to access all Sz sectors, breaking Sz symmetry. In addition, a projection scheme is used to restore global Sz, which gives rise to the Sz spin-projected generalized cluster mean-field (SzGcMF). Both of these extensions contribute to accounting for inter-cluster correlations. We benchmark these methods on the 1D, quasi-2D, and 2D J1 - J2 and XXZ Heisenberg models. Our results indicate that the new methods (GcMF and SzGcMF) provide a qualitative and semi-quantitative description of the Heisenberg lattices in the regimes considered, suggesting them as useful references for further inter-cluster correlations, which are discussed in this work.

Robust formulation of Wick's theorem for computing matrix elements between Hartree-Fock-Bogoliubov wavefunctions.

Numerical difficulties associated with computing matrix elements of operators between Hartree-Fock-Bogoliubov (HFB) wavefunctions have plagued the development of HFB-based many-body theories for decades. The problem arises from divisions by zero in the standard formulation of the nonorthogonal Wick's theorem in the limit of vanishing HFB overlap. In this Communication, we present a robust formulation of Wick's theorem that stays well-behaved regardless of whether the HFB states are orthogonal or not. This new formulation ensures cancellation between the zeros of the overlap and the poles of the Pfaffian, which appears naturally in fermionic systems. Our formula explicitly eliminates self-interaction, which otherwise causes additional numerical challenges. A computationally efficient version of our formalism enables robust symmetry-projected HFB calculations with the same computational cost as mean-field theories. Moreover, we avoid potentially diverging normalization factors by introducing a robust normalization procedure. The resulting formalism treats even and odd number of particles on equal footing and reduces to Hartree-Fock as a natural limit. As proof of concept, we present a numerically stable and accurate solution to a Jordan-Wigner-transformed Hamiltonian, whose singularities motivated the present work. Our robust formulation of Wick's theorem is a most promising development for methods using quasiparticle vacuum states.

State Preparation of Antisymmetrized Geminal Power on a Quantum Computer without Number Projection.

The antisymmetrized geminal power (AGP) is equivalent to the number projected Bardeen-Cooper-Schrieffer (PBCS) wave function. It is also an elementary symmetric polynomial (ESP) state. We generalize previous research on deterministically implementing the Dicke state to a state preparation algorithm for an ESP state, or equivalently AGP, on a quantum computer. Our method is deterministic and has polynomial cost, and it does not rely on number symmetry breaking and restoration. We also show that our circuit is equivalent to a disentangled unitary paired coupled cluster operator and a layer of unitary Jastrow operator acting on a single Slater determinant. The method presented herein highlights the ability of disentangled unitary coupled cluster to capture nontrivial entanglement properties that are hardly accessible with traditional Hartree-Fock based electronic structure methods.

Thermofield Theory for Finite-Temperature Electronic Structure.

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. They also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

Strong-weak duality via Jordan-Wigner transformation: Using fermionic methods for strongly correlated su(2) spin systems.

The Jordan-Wigner transformation establishes a duality between su(2) and fermionic algebras. We present qualitative arguments and numerical evidence that when mapping spins to fermions, the transformation makes strong correlation weaker, as demonstrated by the Hartree-Fock approximation to the transformed Hamiltonian. This result can be rationalized in terms of rank reduction of spin shift terms when transformed to fermions. Conversely, the mapping of fermions to qubits makes strong correlation stronger, complicating its solution when one uses qubit-based correlators. The presence of string operators poses challenges to the implementation of quantum chemistry methods on classical computers, but these can be dealt with using established techniques of low computational cost. Our proof of principle results for XXZ and J1-J2 Heisenberg (in 1D and 2D) indicates that the JW transformed fermionic Hamiltonian has reduced complexity in key regions of their phase diagrams and provides a better starting point for addressing challenging spin problems.

Coupled Cluster and Perturbation Theories Based on a Cluster Mean-Field Reference Applied to Strongly Correlated Spin Systems.

We introduce perturbation and coupled-cluster theories based on a cluster mean-field reference for describing the ground state of strongly correlated spin systems. In cluster mean-field, the ground state wave function is written as a simple tensor product of optimized cluster states. The cluster language and the mean-field nature of the ansatz allow for a straightforward improvement which uses perturbation theory and coupled-cluster to account for intercluster correlations. We present benchmark calculations on the 1D chain and 2D square J1-J2 Heisenberg model, using cluster mean-field, perturbation theory, and coupled-cluster. We also present an extrapolation scheme that allows us to compute thermodynamic limit energies accurately. Our results indicate that, with sufficiently large clusters, the correlated methods (cPT2, cPT4, and cCCSD) can provide a relatively accurate description of the Heisenberg model in the regimes considered, which suggests that the methods presented can be used for other strongly correlated systems. Some ways to improve upon the methods presented in this work are discussed.

A power series approximation in symmetry projected coupled cluster theory.

Projected Hartree-Fock theory provides an accurate description of many kinds of strong correlations but does not properly describe weakly correlated systems. On the other hand, single-reference methods, such as configuration interaction or coupled cluster theory, can handle weakly correlated problems but cannot properly account for strong correlations. Ideally, we would like to combine these techniques in a symmetry-projected coupled cluster approach, but this is far from straightforward. In this work, we provide an alternative formulation to identify the so-called disentangled cluster operators, which arise when we combine these two methodological strands. Our formulation shows promising results for model systems and small molecules.

Construction of linearly independent non-orthogonal AGP states.

We show how to construct a linearly independent set of antisymmetrized geminal power (AGP) states, which allows us to rewrite our recently introduced geminal replacement models as linear combinations of non-orthogonal AGPs. This greatly simplifies the evaluation of matrix elements and permits us to introduce an AGP-based selective configuration interaction method, which can reach arbitrary excitation levels relative to a reference AGP, balancing accuracy and cost as we see fit.

Exploring non-linear correlators on AGP.

Single-reference methods such as Hartree-Fock-based coupled cluster theory are well known for their accuracy and efficiency for weakly correlated systems. For strongly correlated systems, more sophisticated methods are needed. Recent studies have revealed the potential of the antisymmetrized geminal power (AGP) as an excellent initial reference for the strong correlation problem. While these studies improved on AGP by linear correlators, we explore some non-linear exponential Ansätze in this paper. We investigate two approaches in particular. Similar to Wahlen-Strothman et al. [Phys. Rev. B 91, 041114(R) (2015)], we show that the similarity transformed Hamiltonian with a Hilbert-space Jastrow operator is summable to all orders and can be solved over AGP by projecting the Schrödinger equation. The second approach is based on approximating the unitary pair-hopper Ansatz recently proposed for application on a quantum computer. We report benchmark numerical calculations against the ground state of the pairing Hamiltonian for both of these approaches.

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems.

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. The resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

Correlating the antisymmetrized geminal power wave function.

Strong pairing correlations are responsible for superconductivity and off-diagonal long-range order in the two-particle density matrix. The antisymmetrized geminal power wave function was championed many years ago as the simplest model that can provide a reasonable qualitative description for these correlations without breaking number symmetry. The fact remains, however, that the antisymmetrized geminal power is not generally quantitatively accurate in all correlation regimes. In this work, we discuss how we might use this wave function as a reference state for a more sophisticated correlation technique such as configuration interaction, coupled cluster theory, or the random phase approximation.

Geminal Replacement Models Based on AGP.

The antisymmetrized geminal power (AGP) wave function has a long history and is known by different names in various chemical and physical problems. There has been recent interest in using AGP as a starting point for strongly correlated electrons. Here, we show that in a seniority-conserving regime, different AGP-based correlator representations based on generators of the algebra, killing operators, and geminal replacement operators are all equivalent. We implement one representation that uses number operators as correlators and has linearly independent curvilinear metrics to distinguish the regions of Hilbert space. This correlation method called J-CI provides excellent accuracy in energies when applied to the pairing Hamiltonian.

Efficient evaluation of AGP reduced density matrices.

We propose and implement an algorithm to calculate the norm and reduced density matrices (RDMs) of the antisymmetrized geminal power of any rank with polynomial cost. Our method scales quadratically per element of the RDMs. Numerical tests indicate that our method is very fast and capable of treating systems with a few thousand orbitals and hundreds of electrons reliably in double-precision. In addition, we present reconstruction formulas that allow one to decompose higher order RDMs in terms of linear combinations of lower order ones and geminal coefficients, thereby reducing the computational cost significantly.

Thermofield Theory for Finite-Temperature Coupled Cluster.

We present a coupled cluster and linear response theory to compute properties of many-electron systems at nonzero temperatures. For this purpose, we make use of the thermofield dynamics, which allows for a compact wave function representation of the thermal density matrix, and extend our recently developed framework ( J. Chem. Phys. 2019 , 150 , 154109 , DOI: 10.1063/1.5089560 ) to parametrize the so-called thermal state using an exponential ansatz with cluster operators that create thermal quasiparticle excitations on a mean-field reference. As benchmark examples, we apply this method to both model (one-dimensional Hubbard and Pairing) and ab initio (atomic Beryllium and molecular Hydrogen) systems, while comparing with exact results.

Thermofield theory for finite-temperature quantum chemistry.

Thermofield dynamics has proven to be a very useful theory in high-energy physics, particularly since it permits the treatment of both time- and temperature-dependence on an equal footing. We here show that it also has an excellent potential for studying thermal properties of electronic systems in physics and chemistry. We describe a general framework for constructing finite temperature correlated wave function methods typical of ground state methods. We then introduce two distinct approaches to the resulting imaginary time Schrödinger equation, which we refer to as fixed-reference and covariant methods. As an example, we derive the two corresponding versions of thermal configuration interaction theory and apply them to the Hubbard model, while comparing with exact benchmark results.

Polynomial-product states: A symmetry-projection-based factorization of the full coupled cluster wavefunction in terms of polynomials of double excitations.

Our goal is to remedy the failure of symmetry-adapted coupled-cluster theory in the presence of strong correlation. Previous work along these lines has taken us from a diagram-level analysis of the coupled-cluster equations to an understanding of the collective modes which can occur in various channels of the coupled-cluster equations to the exploration of non-exponential wavefunctions in efforts to combine coupled-cluster theory with symmetry projection. In this manuscript, we extend these efforts by introducing a new, polynomial product wavefunction ansatz that incorporates information from symmetry projection into standard coupled-cluster theory in a way that attempts to mitigate the effects of the lack of size extensivity and size consistency characteristic of symmetry-projected methods. We describe the new approach in detail within the context of our previous efforts, explore some illustrative calculations, and consider one route for reducing the computational cost of the new method.